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Abstract 

A method for calculating multi-portfolio time consistent multivariate risk measures 
in discrete time is presented. Market models for d assets with transaction costs or 
illiquidity and possible trading constraints are considered on a finite probability space. 

The set of capital requirements at each time and state is calculated recursively back¬ 
wards in time along the event tree. We motivate why the proposed procedure can be 
seen as a set-valued Bellman’s principle, that might be of independent interest within 
the growing field of set optimization. We give conditions under which the backwards 
calculation of the sets reduces to solving a sequence of linear, respectively convex vector 
optimization problems. Numerical examples are given and include superhedging under 
illiquidity, the set-valued entropic risk measure, and the multi-portfolio time consistent 
version of the relaxed worst case risk measure and of the set-valued average value at 
risk. 

Keywords and phrases: dynamic risk measures, transaction costs, set-valued risk 
measures, set optimization, vector optimization, algorithms, dynamic programming, 
Bellman’s principle 
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1 Introduction 

Multivariate risk measures, also called set-valued risk measures, were introduced in a 
static one-period setting by Meddeb, Jouini, Touzi m in 2004, and further studied 
in puna eh to consider risk evaluations of random vectors motivated by market 
models with transaction costs. Dynamic set-valued risk measures were presented in 
@0 0 [2]. The set-valued version of time-consistency, called multi-portfolio time 
consistency, was introduced in [7J and it was shown to be equivalent to a recursive 
form for the risk measure. 
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In the present paper we will show that this recursive form can be seen as a set¬ 
valued version of Bellman’s principle. On one hand, it enables us to calculate the value 
of a risk measure, that is, the set of all risk compensating initial portfolio holdings, 
backwards in time. This is in the spirit of dynamic programming. On the other hand, 
one can show that the principle of optimality holds true: the truncated optimal strategy 
calculated at time t = 0 is still optimal for the optimization problems appearing at any 
later time point t > 0. 

We will provide conditions under which the recursive form is equivalent to a se¬ 
quence of linear vector optimization problems, which can be solved with Benson’s 
algorithm, see e.g. PH EBII21 El- This is the case for most coherent, but also for con¬ 
vex polyhedral risk measures. More generally, we will give conditions under which 
the recursive form is equivalent to a sequence of convex vector optimization problems, 
that can be approximately solved by the algorithms proposed in [21] , Numerical ex¬ 
amples will illustrate the results. We will give a financial interpretation of the optimal 
strategies, which provides an intuition for the nested formulation of the risk measures. 

The organization of this paper is as follows. In section [2] we will give a short 
introduction to dynamic set-valued risk measures and multi-portfolio time consistency 
based on The main results of this paper are in sections [3] and [4j where we 

reformulate the problem in terms of vector optimization. We will discuss the linear 
case in section |4.1| and the case of convex vector optimization problems in section 
In section [5] we give interpretations and the link to Bellman’s principle of optimality. 
In section [6j we study the calculation of market extensions of set-valued dynamic risk 
measures that appear when trading opportunities at the market are taken into account. 
Interpretations of the results will be given. Numerical examples will be studied in 
section [7j including superhedging, the relaxed worst case risk measure, the average 
value at risk, and the entropic risk measure. 


4.2 


2 Set-valued dynamic risk measures 

Consider a discrete time setting T = {0,1,...,T} with finite time horizon T and a 
finite filtered probability space (D, J 7 , (Ft)teri P) with the power set sigma algebra, 
i.e. T = 2 n , and J~t = J~. Denote by L t = L°(0, J 7 ), P; M d ) the linear space of the 
equivalence classes of Tj-measurable functions X : D —> and by L := L T . 

We write L t , = {X € L t : X € P-a.s.} for the closed convex cone of Revalued 

Ju-measurable random vectors with non-negative components. Note that an element 
X € L t has components Xi,...,Xd in L t (M) = L°(f2, J-), P; R). More generally, we 
denote by L t (Df ) those random vectors (or variables) in L t (resp. L t (M)) that take 
P-a.s. values in Dt■ 

As in [16] and discussed in mm, the portfolios in this paper are in ‘physical 
units’ of an asset rather than the value in a fixed numeraire. That is, for a portfolio 
I 6 tj, the values of Xi (for 1 < * < d) are the number of units of asset i in the 
portfolio at time t. 

Let the set of eligible portfolios M C i.e. those portfolios which can be used to 
compensate for the risk of a portfolio, be the same for all times t and states well, Let 
M be a finitely generated linear subspace of R^ with M + := M n ^ {0}. Denote 
Mf ■= L t (M) and M* + = L t (M + ). Of particular interest, especially when dealing with 
the market extension discussed in section [6j is the case where all assets are eligible, 
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i.e., M = M d . 

A conditional risk measure is a function which maps a d-dimensional random vari¬ 
able into := {D C M t : D = D + M t+ }. a subset of the power set of 

Mt . That is, conditional risk measures map into collections of random vectors in Mt. 
Conditional risk measures were defined in [7] as follows. 

A function Rt : L —>• V{Mt ; Mt,+) is a conditional risk measure at time t if it is 

1. Alt- translative: Vmt G Mt '■ Rt{X + mt) = Rt.(X) — mt ; 

2. L + -monotone: Y — X G L + =>■ Rt(Y) D Rt{X ); 

3. finite at zero: 0 / -Rt(O) / M t . 

A conditional risk measure is called normalized if 

= tftpo + ^i(o) 

holds for every I 6 Lj. It is local if for every A G A and every A G Ju it holds 

lyl-RtPO = 

where In is the indicator function. A conditional risk measure at time t is called 
conditionally convex if for all X, Y G L, and for all A G L t ([ 0,1]) 

R t (XX + (1 - X)Y) D A Rt(X) + (1 - A )R t (Y), 

it is conditionally positive homogeneous if for all X G L and for all A G L t (R ++ ) 

Rt{ A A) = XRt(X), 

and it is conditionally coherent if it is conditionally convex and conditionally posi¬ 
tive homogeneous. A conditional risk measure Rt is closed if graph(i?j) := {(A, u) G 
L x Mt : u G Rt( A)} is closed in the product topology. It is called (conditionally) 
convex upper continuous if 

R~ 1 (D) := {A G L : R t {X) nD/0} 

is closed for any closed (conditionally) convex set D G Q{Mp, Mt-) := {D C Alt '■ D = 
cl co (D + M t -)}. 

A dynamic risk measure is a sequence of conditional risk measures (Rt)teT- It 
is said to have one of the above properties if Rt has this property for every t G T. The 
acceptance set at time t associated with a conditional risk measure Rt is defined by 

A t = {A G L : 0 G R t { A)} . 

Let us now discuss the economic interpretation of a dynamic risk measure {Rt)t£ t- 
At time t = 0, one would choose a portfolio uq G Rq(X), usually to be as small as 
possible, that is, a weakly minimal element of the set to be put aside to make 

A acceptable at time T according to the acceptance set Aq. As time progresses to 
t = 1 and new information become available, one would update this risk compensating 
portfolio to keep the overall position acceptable (now according to A \) by choosing a 
portfolio u\ G — uq + R.\ (A), again usually a weakly minimal element of this set. This 
could mean injecting more capital/assets or extracting them. Now in total uq + u\ G 
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-Ri(X) has been put aside to compensate the risk of X. This procedure continues until 
time T. Thus, at any time 1 £ T one has put aside a portfolio that compensates for 
the risk of X according to the time t risk measure Rt- 

To ensure that updating the risk evaluation is done in a time consistent way, the 
concept of time consistency for scalar risk measures was extended to the set-valued 
framework in [7|| and is called multi-portfolio time consistency. 

Definition 2.1. A dynamic risk measure (Rt)teJ called multi-portfolio time consis¬ 
tent if for all times t £ T\{T}, all portfolios X £ L, and all sets Y C L the implication 

Rt+i(X) C \J Rt+i(Y) => Rt.(X) C \J R t (Y) (2.1) 

Ye Y Ye Y 


is satisfied. 

In |3 theorem 3.4], it was shown that for a normalized dynamic risk measure 
(Rt)t£ Ti with R t : L —> V(M t ; M t:+ ) for all times t, (Rt)teJ being multi-portfolio time 
consistent is equivalent to (Rt)te t being recursive, that is for every time t £ T\{T} 

Rt(X)= J Rt,t+i(-Z) =: R t , t+ i(-R t+ 1 (X)), (2.2) 

ZcR t+ 1 {X) 

where Rt,t+i ■ Mt+i ~> M ij+ ) denotes the stepped risk measure Rt,t+i = Rt\M t+1 , 

that is the restriction of Rt to Mt+ 1 - 

Given an arbitrary dynamic risk measure (Rt)teJ on L, one can compose the one- 
stepped risk measures Rt,t+i backwards in time to obtain a multi-portfolio time con¬ 
sistent risk measure (Rt)t£ T as follows: For all X £ L define 

Rt(X) = R t {X), (2.3) 

Vt £ {0,1,..., T — 1} : Rt(X)= J R tjt+ 1 (-Z). (2.4) 

ZeRt+RX) 


Then, (Rt)te.J is multi-portfolio time consistent and satisfies the properties of M t - 
translativity and monotonicity of a dynamic risk measures, but may fail to be finite 
at zero. Additionally, if (Rt)te t is (conditionally) convex, (conditionally) coherent, or 
(conditionally) convex upper continuous and (conditionally) convex, then (Rt)te t has 
the same property, see proposition 3.11 in fTj and proposition 5.1 in [9j. 

Example 2.2. The set-valued average value at risk was introduced and computed in 
in the static setting. In w the definition was extended to the dynamic framework as 
follows. For parameters A* £ L t , where 0 < A* < 1 and X £ L define 

AV@Rt(X) = {diag (Ai) -1 E [Z\ X t ] ~ z : Z £ L + , X + Z - z £ L + , z £ L t } n M t . 

(AV@Rt)t£j is a normalized closed conditionally coherent dynamic risk measure. The 

multi-portfolio time consistent version AV@R t was discussed in w, with the dual 
representation deduced for the case M = M rf . Since we are only considering a finite 
probability space, we can immediately conclude (as in m that the dynamic average 
value at risk is a polyhedral risk measure. 
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3 A set-valued Bellman’s principle 


We now want to answer the question if it is possible to use the nested formulation (2.2), 
or, more generally, the backward composition (2.3), (2.4) to explicitly calculate the set 
Rt(X), respectively the multi-portfolio time consistent version Rt(X), backwards in 
time. If this is possible it would justify calling this procedure a set-valued Bell¬ 
man’s principle, yielding a dynamic programming method for set-valued functions. 
This would be an interesting insight in itself within the field of set-optimization with 
applications beyond the one considered here. 

Recall that we assumed a finite sample space II with the power set sigma algebra, 
i.e. T = 2 n . We define fit as the set of atoms in T t - For any oj t G fit we denote the 
successor nodes by 

succ(wi) = {ut+i G fl t+ 1 : uj t+ 1 C u t }. 


We use the convention that for an Tt-measurable random variable u, we denote by 
u(ujt) the value of u at node u>t, that is u{uit) ■= u(u) for some uj G w*. Further, we 
denote by Rt(X)[ut\ := : u G Rt(X)} the collection of projections of elements 

of Rt{X) onto 0Jt- Though Rt(X) is a collection of random variables rather than a 
random set, Rj(V)[-] is a random set. 

In order to study a possible calculation of a set of random variables Rt(X) back¬ 
wards in time on a finite event tree, one first needs to check if one can calculate the 
set Rt(X) u)t- wise at each node. That is, we wish to verify that 


u G R t (X) u{u t ) G Rt{X)[u> t \ G fl t . 


This consideration is not a concern when dealing with scalar risk measures, but in the 
set-valued case certain conditions are needed to ensure it. Since we work on a finite 
probability space, for an uj t -wise approach to (2.3), (2.4) one will need {Rt)te T to have 
^-decomposable images, i.e., 


Rt(X) = 1 ARt{X) + l AcRt(X) VX g L VA g J), 


which is satisfied if (Rt)te t has conditionally convex images. Furthermore, for the 
multi-port folio time consistent version to only depend on the possible future (successor) 
states, one will need (Rt)te T to be local. 

Remark 3.1. Assuming (Rt)tej to be conditionally convex implies both, (Rt)te.J be¬ 
ing local (see proposition 2.8 in 0), as well as (Rt)te T being conditionally convex 
(proposition 3.11 in [Jj) and thus having conditionally convex images. 


In order to implement the set optimization problem, we wish to reframe the back¬ 
ward composition (2.3), (2.4) as a vector optimization problem. For this reason we 
require that (Rt)te T has closed images as well. 


Remark 3.2. It is challenging to ensure that (Rt)te T has closed images. Let us give 
three examples, where (Rt)te t is closed (and thus has closed images): a) if the dynamic 
risk measure {Rt)te T is convex upper continuous and convex (see proposition 5.1 in [9]), 
or b) if the dynamic risk measure (Rt)t £T is conditionally convex upper continuous and 
conditionally convex, or c) if the dynamic risk measure (Rt)te T is polyhedral (that 
is if graph(Rf) is a convex polyhedron, i.e. the intersection of finitely many closed 
half-spaces). 
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If (Rt)teT defined in (2.3), (2.4) is conditionally convex, but does not already have 


closed images, one needs to consider its closed-valued version, i.e. 


Vt E {0,1,..., T — 1} : 


R t {X) := cICRtPO) 

Rt{X) := cl [J Rt(-Z) 

Z£R t+ i(X) 


(3.1) 

(3.2) 


for all portfolios X E L. We will show that (Rt)t£j can be approximated for arbitrarily 
small 5 > 0 by (Rt)te T that admits an u^-wise representation. The approximation is 
understood in the following sense. {Rt)teJ is called an approximation of (Rt)teT if 

R t (X) + Sml C R t (X) C R t (X) 

for any time t, for every portfolio X E L, for any approximation tolerance 6 > 0 and any 
m E int(M + ) in the subspace topology, i.e. with the topology tm '.= {A n M : A £ t} 
where r is the topology on M rf . Thus, if (Rt)te T does not have closed images, we will 


use construction (3.1), (3.2) to calculate an approximation (Rt)teJ of (Rt)te t- 


Theorem 3.3. Let (Rt)te t be a conditionally convex dynamic risk measure. Let 


(Rt)t&J denote its multi-portfolio time consistent version as defined in (2.3), (2.4). 


Then, we can calculate an approximation (Rt)te t of (Rt)te T in on ut-wise manner by 

Rt{X) = cl(i?j i (X)[o;7’]), Vt ut £ 12 t ; (3-3) 

R t {X)[u t \ =d\J{Rt,t + i(-Z)[u t ] : Z(uj t+ i) E R t+ i(X)[uj t+l ] Vu t+ i E succ(o;t)} 

(3.4) 


for every state ut E Lit ond all times t E T\{T}. 


Proof. Define (Rt)te T by equations (3.1), ( |3.2| ). Clearly, (Rt)te T has closed images. One 
can show that (Rt)teT is conditionally convex by backward induction. The assertion 
is clearly true for Rt. Now assume Rt+i is conditionally convex and let X, Y E L and 
A E L t ([ 0,1]). Then, 


XR t (X) + (1 - X)R t (Y) = Acl (j R t (—Z x ) + (1 — A) cl |J Rt(-Z y ) 

z x eR t +i(x) z Y eRt+i(Y) 

Cell cl |J XRtf—Zx) + cl |J (1 - \)R t (-Z Y ) I 

\ z x eRt+i(x) z Y eR t+1 (Y) J 

= cl [J [A Rt(-Zx) + (1 - X)R t (-Z y )\ 

z x eR t+1 (x) 

Z Y £Rt+i (Y) 

c cl J Rt(-(XZx + (1 - A )Z Y )) 

z x eR t+ i(A') 

Z Y SRt+i(Y) 

= Cl U i?t(-Z) 

ze\Rt+i(x)+(i-\)R t+1 (Y) 

C cl J iit(-Z) =Rt(XX + (l-X)Y). 

ZeR t+1 (\X+(l-\)Y) 
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Since the probability space is assumed to be finite, fit is by definition the finest partition 
of 17 in Tt and 1 Ut u G 1 Ut Rt(X) if and only if u(u>t) G Rt(X)[ut\. Since (Rt)teJ has 
conditionally convex images and the probability space is finite, it follows that u G R-t(X) 
if and only if u{ut) G Rt(X)[ujt] for every ut G fit- Thus, one can calculate Rt{X) ut- 


wise. Therefore the terminal condition (3.3) holds trivially. Let t G T\{T} and ut G fit, 


using (3.2) and the local property for Rt (see [H proposition 2.8]) it follows that 


l Ut Rt(X) = l wt cl Rt(-Z) 

zeRt+RX) 

= cl (^J \ Ut Rt(}-u t (—Z)) = cl (^J 1 u t Rt{—Z). 

ZGRt+ i(X) Z£l UJt Rt-\-i(X) 


Note that Z G l Ut Rt+i(X) if and only if Z(ut+i) G ^ + i(X)[o;t + i] for every u>t +1 £ 
succ(wt) and Z(uj t +x) = 0 otherwise. But as we only need to consider 1 Wt Z by Rt local, 
the only constraints on Z are imposed by u>t+i G succ(wi). Thus, (3.4) follows. 

Finally we will show that Rt is an approximation of Rt. By definition, Rt(X) := 
cl (Rt(X)) for all X G L, which implies that Rt{X) + SttuI C Rt{X) C Rt(X) for 
any portfolio X G L, for any St > 0, and any m G int(M+). For the proof by induction, 
assume that is an approximation of Rt+i- Then for t <T — 1 we have 


Rt(X)= U Rt{-Z)C cl U R t (-Z) = R t (X) 
zeR t+1 (x) zeRt+RX) 

C cl Rt{~Z + St+iml) = cl R t (-Z) - St+iml 

zeR t+ 1 {x) zsRt+RX) 

- U Rt(-Z) - (e +S t+ x)ml = Rt(X) - (e +6t+i)ml 

zeRt+RX) 


for any e, 8t+\ > 0, and any m G int(M+). The first inclusion on the second line 
follows from the induction hypothesis. The first inclusion on the third line follows 
since cl (Rt(X)) C R t (X) — em, for any e > 0 and all X G L. Denote St := e + > 0. 

Therefore for any time t and any St > 0 we have that Rt{X) + Sf/rril C Rt(X) C 
Rt{X). □ 


In the recursion (3.4), the calculation is dependent on R t:t +i(—Z)[ut] which is by 
locality of Rt equal to Rt,t+i(— l Wt Z’)[wt], i.e. R^t+i(—Z)[u}t] only depends on the part 
of Z that can be attained from state ut- For a local risk measure, we can therefore 
define i2t,t+i(■)[&;*] on Mt+i[u>t\ := L°(U succ(u;i), 2 succ ^ t \P(-|o;t); M) (the equivalence 
class of T) + i-measurable random variables in the eligible space M starting from node 
ut) by R ty t+i{Z)[u t ] := Rt,t+i{Z)[u t } for Z G M t+ 1 [ut] and some Z G M t+ 1 with 
Z{u t+ 1 ) = Z(u t+ 1 ) for every u t +\ G succ(o;t). 


Remark 3.4. If (Rt)ter defined in (2.3), (2.4) does have closed images already, (Rt)tej 


coincides with (Rt)teJ and thus can be calculated in an ut -wise manner by (3.3), 


Observe that 


is a set-valued optimization problem in the complete lattice 


Q(M t ; Mt,+) '■= {D C Mt : D = clco(Z? + Mi ! -|-)}, and (3.4) is a set-valued optimization 
problem in the complete lattice Q{M ; M+). As such, while the initial proposed problem 


(3.2) requires a lattice in sets of random vectors, the wt-wise representation allows us 
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to consider the lattice in sets of real-valued vectors instead. That is, the objective 
function Rt,t +1 at node ut is a set-valued function that is minimized over the constraint 
set Z t := {Z G M t+1 [u t \ ■ Z(uj t+ i) G Rt+i(X)[u t + 1] Vw t+ 1 G succ(wt)}: 


fl t (X)M = cl M inf i? M+1 (-Z)[u; f ]. 

zez t 


(3.5) 


Recall from [19] that the (lattice) infimum over a function / : [u;t] —>• G(M] M + ) is 

given by infxeA” f(X) = clco|J YgA , f(X ) for X C M i+ i[wt], where the convex hull in 
front of the union can be dropped here as Rt(X)[ut\ is convex by theorem 
the idea of 


one can transform the set-valued problem (3.5) into a linear vector 


3.3 


Using 


optimization problem that can be solved e.g. by Benson’s algorithm if the risk measure 
(Rt)teJ is polyhedral, i.e., the graph of R t is a convex polyhedron. More generally, if 
(Rt)teJ is the upper image of a convex vector optimization problem, one can transform 


the set-valued problem (3.5) into a convex vector optimization problem that can be 
approximately solved by the algorithms discussed in [21] . In both cases, one uses that 


the value of the set-optimization problem (3.5) can be written as the value of a vector 
optimization problem 


Rt(X)[ut\= inf *t{Z,Y) (3.6) 

(z,Y)ez t 

for the linear vector-valued function = {T}, feasible set 

Zt = {(Z,Y) G M t+1 [u t ] x M : 

Y G R t ,t+i{~Z)[ujt],Z(u;t + i) G R t+ i(X)[u) t +i] Vut+i G succ(w t )} 

and ordering cone M + . Let = {d> 4 (Z) : Z G Z t } denote the image of the feasible 

set. The set cl(<h f (.Z)) + M + ) is called the upper image of the vector optimization 
problem (3.6). We now discuss the constraints by looking at two cases: the polyhedral 
case and the convex case. 


4 Computational procedure 

4.1 Linear vector optimization and polyhedral risk mea¬ 
sures 


Recall that a risk measure Rt is polyhedral if its graph is a convex polyhedron, i.e. 
the intersection of finitely many closed half-spaces. It is also equivalent to Rt having 


a polyhedral acceptance set. For a polyhedral risk measure (Rt.)te T> problem (3.6) is a 
linear vector optimization problem. 

Proposition 4.1. If the dynamic risk measure (Rt)tej is conditionally convex and 
polyhedral, its multi-portfolio time consistent version (Rt)te T? defined in (2.3), (2.4), 


can be calculated uit-wise, where in each node ojt G Q t , t G T\{T}, the linear vector 


optimization problem (3.6) has to be solved. 


Proof. The wt-wise representation follows from theorem |3.3[ Now, let us show that 


problem (3.6) is a linear vector optimization problem. By (Rt)t .£T polyhedral and since 












Rt,t+i maps into Q{M\ M+), -Rt,t+i(—^)[wt] is the upper image of a linear vector 
optimization problem (see remark 5.1 in [T2J), thus 


R t ,t+i{-Z)[w t ] = { Pt(x ) + M + : B t x > b t } 


(4.1) 


for a vector x = (Z, z) that might include some auxiliary variable z, and for ma¬ 
trices Pt and Bt and vectors bt of appropriate dimensions. Then, the constraints 
Y E Rt t t+i{—Z)[uit] can be equivalently written as 


(~Z,Y) E graph R t ,t+i M] 


M J (Y-P t (x))> 0, B t x>b t , 


where the matrix M contains the generating vectors of the positive dual Mf of the 
ordering cone M + . Thus, these constraints are linear. To obtain linearity of the other 
constraints, note that (Rt)teJ polyhedral implies (Rt)te t closed (by definition) and 
(Rt)teJ closed. To see the last implication observe that Rt is polyhedral if and only 
if At is a polyhedron. The acceptance set of Rt is given by At = At,t+ 1 + At+i, see 
[7| corollary 3.14], A t ^+i is a polyhedron since Rt,t+i is polyhedral, by backwards 
recursion we assume At+\ is a polyhedron, and the sum of polyhedra is a polyhedron. 
Therefore At is a polyhedron, which is equivalent to Rt being polyhedral (and thus 
closed as well). Thus, in the polyhedral case, (Rt)te T coincides with (Rt)te T- The 
linearity of the constraints Z E Rt,+i (X) follow by induction. The constraints from the 
terminal condition Z(ojt) E Rt(X)[ujt\ = Rt(X)[ut\ are linear for all ujt E I4t by Rt 
polyhedral. Thus, let us assume the constraints Z{ujt+ 1 ) E .Rt+iPOI^t+i] are linear 
for all (jJt+i E succ(wi) for a given node ut E 14*, then we need to show that Rt(X)[ujt] 


is polyhedral. Since Q is assumed to be finite, problem (3.6) is clearly a linear vector 


optimization problem. By (3.5), (3.6) and Rt : t+i{—Z)[ut\ mapping into Q{M\ M+), we 
have Rt(X)[ut] = $>t(Zt) + M+, which is for finite 14 closed and polyhedral. Thus, 
Rt(X)[idt] is the upper image of the linear vector optimization problem (3.6). □ 


Thus, if one assumes (Rt)te t to be a conditionally convex and polyhedral risk 
measure, Rt(X)[uj t \ is for every t and ut the upper image of a linear vector optimization 
problem, which is polyhedral and can be calculated by Benson’s algorithm (see jT2]). 
Then, the set Rt(X) can be calculated backwards in time, by solving at each node a 
linear vector optimization problem. Note that at time t, the problems for each ut E 14 1 
can be calculated in parallel instead of sequentially which reduces computational time. 
Several examples will be discussed in section [7| 

Benson’s algorithm is an appropriate tool to solve the vector optimization problem 


(3.6) as it takes advantage of the fact that the dimension dim(M) of the image space is 
usually significantly smaller than the dimension d x | succ(wt)| +d+ \z\ of the pre-image 
space, where |succ(wt)| denotes the number of successor nodes of uj t and \z\ denotes 


the dimension of the auxiliary variables in (4.1) 


In practice (especially if M is higher dimensional), when the number of vertices 
of the set R t (X)[ut] is very high, one would calculate an approximation of R t (X)[ujt\ 
having fewer vertices, see remark 4.10 in [ L 12|. Then, for the backward recursion, one 
would need to know how the approximation errors accumulate over time. This will be 


discussed in propositions 4.4 and 4.7 below in a more general framework. 


Example 4.2. The relaxed worst case risk measure was introduced in the static frame¬ 
work in example 5.2 of EM- The idea behind the relaxed worst case risk measure is to 
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modify the worst case risk measure so that portfolios with “small” negative components 
can still be acceptable. 

In this paper we consider the dynamic extension of such a risk measure. By defi¬ 
nition it is a polyhedral and conditionally convex, but not conditionally coherent, risk 
measure; the acceptance set at time t is given by 


A 


RWC 

t 


— (—£ + L _,) n L(G ) 


for some level e E Ml and some finitely generated convex cone G D and G 7 ^ ! rf . 
Note that if G = ore = 0 then the relaxed worst case risk measure is equivalent to 
the worst case risk measure. 

Let (Rt)teJ be the relaxed wor st case risk measure with polyhedral acceptance set 
A c . Then, by proposition 4 ■ 1 one can calculate its multi-portfolio time consistent 
version (Rt)te t oot-wise, where in each node ut. E Lit, t E T\{T}, the linear vector 


optimization problem (3.6) has to be solved. It should be noted that in general Rt 7 ^ Rt, 
i.e. the relaxed worst case risk measure is not multi-portfolio time consistent. 


4.2 Convex vector optimization and conditionally convex 
risk measures 

As the upper image of a convex vector optimization problem can only be calculated by 
a polyhedral approximation yielding an inner as well as an outer approximation with 
respect to some error level e (see e.g. mm), we introduce approximations of sets, 
respectively functions. Throughout, we fix a parameter m E intM + . 

Definition 4.3. Given a set S E V(M; M + ) and an error level e > 0, we call a set 
S e E V(M; M + ) an approximation of S, if 

S e + emCSC S £ . 


Given a set-valued function F \ L —t V(M t ; M t + ) and an error level e > 0, we call the 
function F e : L —^ V(M t ; M tt+ ) an approximation of F if 

F e (X) + eml C F(X) C F e (X) for every X E L. 


In the convex case one can in general only approximately calculate the constraint 
set Rt+i(X) in the backward recursion (3.4), respectively (3.5). Let us study the 
robustness of the set-optimization problem (3.5) to this perturbation of the constraints. 


Proposition 4.4. Lete> 0. Let R% +1 (X)[u}t+i\ be an e-approximation of R t +i(X)[ut+i] 
for each ut+i E succ(wt) ; then Rl(X)[ujt\ defined by 


Rl{X)[u t ] := cl J R t , t+1 (-Z)[uj t }, (4.2) 

ze^f 


with 

Z £ t := {Z E M t+1 [u t } : Z(u) t +i) € R% +1 (X)[u t+ 1 ] \/u t+ i £ succ(wi)}, 


is an e-approximation of Rt(X)[ut] defined in (3.4). 
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Proof. The assumption implies Zf C Z\ C Z t — eml. This together with 
transitivity of Rt yields 

Rt(X)[u; t ]= cl \J Rt,t+i(-Z)[u t ]Cd J R t>t+1 (-Z)[u t \ = R e t (X)[u Jt } 

Z£Z t zez; 

C cl R t ,t+i(-Z)[u t ] = cl IJ R t ,t+i(-Z)[u t ] - em = R t (X)[u t } - em. 

Z-\-ETYl\(zZt ZdZf 


(3.5) and 


Thus, Rt(X)[u>t\ + em C R t {X)[oj t ) C R e t (X)[ujt\, i.e. Rf(X)[ut] is a e-approximation 
of Rt(X)[u t ]. □ 


Next, we discuss under which conditions problem (4.2) is a convex vector optimiza¬ 


tion problem and which additional assumptions are necessary to apply the algorithm 
proposed in [21j to calculate a polyhedral approximation of the upper image of this 
convex vector optimization problem. 


Assumption 4.5. a) Let the objective function in ( 4.2 ) be of the form R t ,t+i(Z)[u t \ 
{4q(Y, z)+M .|_ : gt(Z, z) < 0} for an M + -convex vector function Tt, a component¬ 
wise convex vector function gt and a vector z, all of appropriate and finite dimen¬ 


sions. 


b) Let the function 4-h in a) be continuous, and let the feasible set X := {(Z,z) 
gt(Z , z) < 0} satisfy int L/0. 


Assumption 4.5 a) means that the closure of Rtj+i(Z)[oj t ] is itself the upper image 


of a convex vector optimization problem. 


Proposition 4.6. Let the objective function Rt,t+i(Z)[ujt\ in (4.2) satisfy assump¬ 
tion 4-5 a) and let R e t+l (X)[u>t+ 1] be a polyhedron for each u>t+ 1 E succ^). Then, 
problem (14.21) is a convex vector optimization problem. 


Proof. Similar to (3.6), the set-valued problem (4.2) can be written as a vector opti¬ 
mization problem 


inf Q t (Z,Y), 
(z,Y)ez; 


by setting <&t(Z, Y ) = [Y } (which is a linear vector function), defining the feasible set as 
Z\ = {(Z,Y) E M t+1 [u t \ X M : Y E R t , t+1 (-Z)[uj t ], Z(cj t+1 ) E Rt +1 (X)[u t+1 ] Vw t+ 1 E 
succ(a;t)} and using M + as the ordering cone. The constraints Z(ut+i) E R\ +l (X)[ujt+ 1 ] 
are by assumption linear. Under assumption 4.5 a), the constraints Y E Rt t t+i(~Z)[u>t] 
can be equivalently written as 


m 


(Y, -Z) E graph R t>t+ i[ujt\ 


M T (M-Z, Z ) - Y) < 0, g t (-Z, z) < 0, 


where the matrix M contains the generating vectors of A iff. Z, z) — Y) is a 

component-wise convex vector function since is a M + -convex vector function. Thus, 
these are convex constraints and (4.2) is a convex vector optimization problem. □ 


The additional assumptions 4.5 b) are necessary to ensure that problem (4.2) can 
be (approximately) solved by the algorithms presented in ([21]. In detail 


sumptions |4.5| and if the feasible set X := {(Z,z) 


under as- 
gt(Z,z ) < 0} is compact, [21] 
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theorems 4.9 and 4.14] state that the algorithms in [21] provide an approximation of 
the upper image of (4.2), i.e. a polyhedral approximation of Rt(X), if they termi¬ 


nate. However, the compactness assumption is typically not satisfied in the setting of 
risk measures. In that case EH remark 3 in section 4.3] shows that the algorithms 


presented in m still return an approximation of the upper image of ( |4.2[ ) as long as 
all the scalar optimization problems within the algorithm can be solved and the algo¬ 
rithm terminates. In the example of the set-valued entropic risk measure considered 


in section 4.9, this will indeed be the case. 


Since in general problem (4.2) can only be solved approximately (e.g. by the algo¬ 
rithms in 


), one also need to study how the approximation errors made at different 
time points accumulate over time. 


Proposition 4.7. Let e, 7 > 0. If Rf’ 1 (X)[uot\ is a 7 - approximation of Rl(X)[ujt\ 
defined in (4.2), then Rf 1 (X)[ujt\ is an (e + -approximation of Rt(X)[ujt] defined in 

dm 


Proof. 


Rf' 1 (X)[u)t\ being a 7 -approximation of Rl(X)[ojt] means 


R t '\X)[u t ] + 7™ C R\{X)[u t ] C R^{X)[u t ], 


Proposition 4.4 shows that Rl(X)[oj t ] is an e-approximation of R t (X)[ut], i.e. 


R e t (X)[u; t } + em C R t {X)[u t ] C R e t (X)[u t ]. 
Both chains of inclusions yield 

R e S(X)[u t ] + (e + 7 )m C Rt(X)[u t ] C Rf 1 (X)[u t \. 


□ 


We are now ready to prove the main result of this section. Recall that the aim was 
to (approximately) calculate the multi-portfolio time consistent risk measure (Rt)teJ 
backwards in time in the spirit of a set-valued Bellman’s principle. We will see that 
(Rt)teJ can be obtained by solving at each node backwards in time a convex vector 
optimization problem. In practice, these problems can only be approximately solved. 
But we are able to determine the overall approximation error, when the approximation 
error at each node is chosen to be e > 0. One could of course also vary this error level 
at different nodes or different time points and obtain corresponding results. 


Proposition 4.8. Let (Rt)teJ be a conditionally convex dynamic risk measure sat¬ 
isfying assumption |^.5[ Let e > 0. Then for any time t and given X E L, we can 
find a [(T — t + l)e + 5}-approximation of the multi-portfolio time consistent version 
(Rt(X))t£T defined in (2.3), (2.4), by calculating backwards in time at each nodeut E Lit 


an e-approximation of the upper image of the convex vector optimization problem (4.2) 
Here 6 > 0 can be chosen arbitrarily small. 


Proof. Ass ump tion 4.5 and the local property of (Rt)te t imply that all the assumptions 


of theorem 3.3 are satisfied, thus a ^-approximation (Rt)te t of (Rt)te t can be calculated 


wt-wise for arbitrarily small 6 > 0 by (3.4), (3.3). 

For t = T one obtains an e-approximation Rtp(X) of Rt(X) = cl(ii/r(A")) = 
c!(Rt(X)) by calculating an e-approximation of the upper image of the convex vector 
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optimization problem Rt(X)[ljt] = {^t(X,z) + M + : gx(X,z) < 0} (see assump¬ 
tion 4.5 a)) at each node ujt E fiy- R e T (X) is by construction polyhedral (see the 
algorithms in [2T]) and is the input for problem (4.2) at time t = T — 1 , which is 
by proposition |4.6| then a convex vector optimization problem. Its solution would by 
proposition 


4.4 


yield an e-approximation R f T _\ (X) of Rt~i(X), but one can in general 
only calculate an e-solution. This e-solution yields an e-approximation of Rf,_ 1 (X), 
which is by proposition 4.7 a 2e-approximation of Rt~i(X). 

Going backwards like this yields for any t a (T — t + l)e-approximation of Rt(X), 
which is by theorem |3.3| and the logic of adding up approximation errors as in proposi¬ 


tion 4.7 a (T — t + l)e + ^-approximation of the multi-portfolio time consistent version 
Rt(X) for arbitrarily small 5 > 0. □ 

Example 4.9. The set-valued, entropic risk measure was studied in m in a single 
period static framework. The dynamic version was discussed in m ■ As in the scalar 
case, the entropic risk measure is intimately related to the exponential utility function. 
Consider risk aversion parameters X E T t (M^. + ), Ct E Q{L t \L t+ ) with 0 £ Q and 
Ct n = 0. The dynamic entropic risk measure is defined by 


R e t nt (X- X, C t ) := {u E Mt : E [u t (X + u)\ Ft] E C t } 

, X , / N N T 


(4.3) 


for every X E L where ut{x) = (ut,i(x\ ),..., ut,d{xd)) for any x E Mr and ut,i(z ) = 


l—e~ 


for zgR and i = 1 ,..., d. 

An approximate calculation of the static entropic risk measure was shown in m via 
solving a convex vector optimization problem. With the method presented in section [3] 
we are able to compute an approximation (Rf nt )t£j of the multi-portfolio time consistent 
version t by backward composition for a general space of eligible portfolios Alt, 

(stochastic) risk aversion parameters X, and polyhedral parameters Ct- It was proven 
in m that the entropic risk measure is c.u.c. and multi-portfolio time consistent in the 
case that M = M. d , constant A E K+_(_, and Ct = L tj +, *-e. Rf nt = R*( nt = Rf nt . 

From the definition of the entropic risk measure with Ct polyhedral, it is clear that 


Rt satisfies assumption 4-5 Thus, by proposition 4-8 one can calculate an approxi 


mation of the multi-portfolio time consistent version (Rf nt (X))t£i by calculating back¬ 
wards in time at each node ut E lit an approximation of the upper image of the convex 


vector optimization problem (4.2) using the algorithms presented in M- 


5 Interpretation and relation to Bellman’s prin¬ 
ciple 


The risk measure (Rt)teT, while constructed backwards in time, has a nice financial 
interpretation involving portfolio injections made as time progresses, that is an inter¬ 
pretation forwards in time: For every choice of a risk compensating portfolio holding 
Zq E Rq(X) at time t = 0, there exists, by equation (2.4), a sequence of portfolio 


holdings (2u)t£T\{o} such that 


and 


Z t e R t {X) 


Zt-\ € Rt.-i{—Zt 


(5.1) 


(5.2) 
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Inclusion (5.1) means Zq, Z±,Zt are the risk compensating portfolio holdings at 


times 0,1, An intuitive interpretation of (5.2) can be obtained by the following 

reformulation. Defining the portfolio injections (respectively withdrawals - if negative) 
( ut)ter that are needed to update the risk compensating portfolio holdings by ut = 
Zt — Zt- i (with uq = Zo), the two conditions on (Zt)t<=T can be rewritten in terms of 
( ut)teJ as follows 

Ut £ Rt 



for every time t , and 


0 £ -Rt,t+i(— ut+i), 


(5.3) 


for t £ T\{T}. Inclusion (5.3) means the risk of the portfolio injection needed at time 
t + 1 is acceptable at time t with respect to the one-period risk measure Rt,t+ 1 - This 
gives the main interpretation of the backward composition of (i?j)t e -j. At each one- 
period step the original measure (Rt)t£T is used, but it is used in a time consistent way 
in the sense of Bellman. 

One can observe Bellman’s principle of optimality: The at t truncated optimal 
solution ( Z s )J =t obtained at time 0 from (2.4) and a given Zo £ Rq(X) is still optimal 
at any later time point t £ T. To see that, note that for the risk compensating 
portfolio holding Zt £ Rt(X), ( Z s )J =t satisfies the conditions Z s £ R S (X) and Z s _i £ 


R s -i(-Z s ), s £ {t ,..., T} from (2.4). 


Let us now explain on how to compute (Rt(X)) t £T and how to obtain for a given 
Zq £ Ro(X) at time 1 = 0 a sequence (Z t )t^ t\{o} of risk compensating portfolio hold¬ 
ings on the realizing path. (Rt)teT can be calculated with the approach discussed in 


sections 4.1 and 4.2 Benson’s algorithm also calculates a solution of the linear vector 
optimization problems in the sense of definition 2.20 in m, respectively, an e-solution 
in the sense of definition 3.3 in [21] for a convex vector optimization problem. These 
finite solution sets are then used to calculate the sequence of risk compensating port¬ 
folio holdings (Z t )teT, respectively the injection/withdrawal strategy (ut)teT, forwards 
in time on the realizing path by solving an additional linear program, specified in the 
following, at each point in time. Let A)^] = {{Z l t+1 [ujt\, Yt[ u t]) '■ i = 1,, n, n £ 


N} C M t +i[u)t\ X M be the (e-)solution set to the vector optimization problem (3.6). 
Let us first explain the method in case of a linear vector optimization problem: For 
any Zo in the risk measure Rq(X), there exists a convex combination of elements of the 
solution on the efficient frontier (the collection of nondominated vectors with ordering 
M + ) such that Zq > A *Yq- This coefficient vector A* £ M" can be found by 

solving any linear optimization problem of the form 


min c T (Yn , ■ 

AeRi v 


, Y 0 n ) A subject to (Tq 1 , • • • , y 0 n ) A < Z 0 , 


1 t A = 1 


(5.4) 


with c £ M^_\{0}. The coefficient vector A* £ K" can then be used to define Zq := 
)T)” =1 A*Yq on the efficient frontier of Rq(X). Notice that Zq = Zq if Zo is already on 
the efficient frontier. Additionally, the next time step full capital requirement is given 
by Z\ := w hich might not be on the efficient frontier of iii(A'). This process 

is repeated through the event tree forwards in time. The choice of cost vector c (or 
alternatively a nonlinear cost function) determines the possible liquidation/withdrawal 
strategy akin to that discussed for the super hedging risk measure in [20] , 
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In the case of a convex vector optimization problem, one can calculate only a poly¬ 
hedral approximation (e.g.with error level e = (T + l)e + <5 as in proposition 4.8) Rq(X) 
of Ro(X). Thus, when choosing the initial capital, one would pick a minimal capital 
from the calculated inner approximation, and not the true set, i.e. uo E Rf(X) + cm. 
Noting that an e-solution of problem (3.6) provides a solution to the linear vector 
optimization problem whose upper image is the inner approximation, the same proce¬ 
dure as in the linear case can be applied, just replacing Rq(X) by its inner polyhedral 
approximation. One obtains an ej-optimal strategy of risk compensating portfolios 
(Z t )t£T, respectively portfolio injections (ut)te t, with e t = (T — t + l)e + 5 when using 
the same error level e > 0 in each iteration step, see proposition |4.8[ 


6 Notes on market extensions 


Market extensions are considered when one is not only interested in putting a ‘capital 
requirement’ u E Rt(X) at time t aside and holding it until time T to make X risk neu¬ 
tral, but in exploiting the trading opportunities at the market to minimize the amount 
of capital needed for risk compensation. For the definition of the market extension be¬ 
low, we will set M = M rf , i.e. we consider the full space of eligible assets. A justification 
for that comes from a mathematical as well as an interpretational aspect, which will be 
detailed in remark 6.4 below. But one can already understand that choice by realizing 
that the role of M comes mainly from a regulatory point of view. A regulator might 
only allow capital requirements to be made in certain currencies for example, and these 
capital requirements are held until time T. But the market extension is more linked to 
internal risk measurement and management as one is exploring trading opportunities 
in possibly all assets, and thus will hold at any time t a portfolio in possibly all assets, 
so there is no need in restricting the capital requirements to be made in certain assets 
only. 

The market extension ( R l t nar )t£T of a dynamic risk measure (Rt)te t is given by 


R r r r (x) ■■= u Mx-k) 

keKt 


for some C L modeling the set of attainable claims. When C L t then it 
immediately follows that 

R™ ar (X) = R t {X) + K t , (6.1) 

Let us give a few examples, all are special cases of the set-valued portfolios introduced 
in j3jj . 


Example 6.1. In a market with proportional transaction costs, trading is modeled by 
a sequence of solvency cones (K t )te t> see m \TflJ . I\ t is a solvency cone at time 
t if it is an Xt.-measurable cone such that for every w G H, K t [u\ is a closed convex 
cone with C Kt[uj\ C M d . Kt is generated by the bid and ask prices between any 
two assets at time t. In a market with proportional transaction costs, one would set 
K t = L t {K t ) (see Q/j. 


Example 6.2. More generally, in markets with illiquidity (convex transaction costs) 
as in f2B/ . trading is modeled by a sequence of convex solvency regions (Kt)t& t- K t 
is a convex solvency region at time t if it is an Tt-tnz&surable set such that for 
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every weO, K t [u\ is a closed convex set with C K t [cv\ C M d . Then, one would set 
Kf = L t (K t ). 

Example 6.3. One could also incorporate trading constraints on the size of trans¬ 
actions by considering convex random sets D t (not necessarily mapping into the closed 
convex upper sets Q(R d ,W^_)) as follows. Given t £ T, let Dt : 0, -> 2 Rd (with 2 m<< denot¬ 
ing the power set ofR d ) be an Tt-measurable function such that D t [uj] is a closed convex 
set and K t [uj\ n D t [uo] / 0 for every u; £ fh Then, one would set K t = L t (K t n D t ). 


We now want to consider the market extended multi-portfolio time consistent ver¬ 
sion of (Rt)teT with respect to C L t . Different possibilities arise and will be briefly 
discussed. First note that the market extension of the multi-portfolio time consistent 
version (Rt)te t > given by ( R t (-) + is not multi-portfolio time consistent, thus 

the solutions do not satisfy Bellman’s principle as detailed in section [5] and there is no 
good economic interpretation. Therefore, only carefully alternating a market exten¬ 
sion step ( 6 . 1 ) and a backward recursion step, that is applying ( 6 . 1 ) in each step of the 
backward recursion will yield the desired results and interpretations. As such we will 
define 


j^rnar 


(X):= |J (R t (-Z) + K t )= |J |J R t (-Z). 


( 6 . 2 ) 


zeR™{ 


(X) 


fceKt zeR™{(x-k) 


This is the multi-portfolio time consistent version of the market extension, i.e., the 
backward composition of (i?™ ar ) te t- The obtained risk measure (R] nar )t £t does indeed 
satisfy Bellman’s principle and we will give the economic interpretation of the solutions 
below. Equation (6.2) shows that the two operations ‘market extension’ and ‘multi¬ 


portfolio time consistent version’ are interchangeable at any single time point t within a 
recursion step under the assumption M = M rf . However, as noted previously, the market 
extension of the multi-portfolio time consistent version, (Rt(-) + Kt)t £ T, is not equal to 
the multi-port folio time consistent version of the market extension, (i? 5 7 iar ) tg -j, defined 
in (6.2). Thus, intrinsic to the definition of (Rf iar )t £t is that both the market extension 


and backward recursion (independent of the order of operations) are computed at a 
time point t + 1 and used as the input for the backward recursion in the next time 
point t. 

In analogy to section [5] for the regulator risk measure, one can obtain a nice finan¬ 
cial interpretation of the market extended composed risk measure (i?™ ar )t e -j involving 
portfolio injections and trades made forwards in time. In addition to the sequence 
of portfolios holdings (Zt)t£j one obtained for the regulator risk measure, there addi¬ 
tionally exists by equation ( |6.2[ ) a sequence of trades (kt)teT such that Z t £ R 7 f Lar {X ), 
kt £ K t, and Z t —kt £ Rt{—Z t + 1 ) for every choice of a risk compensating portfolio hold¬ 
ing Zq £ R™ ar (X) at time t = 0. That means Zq, Z\, .... Zt are the risk compensating 
portfolio holdings before trades at times 0,1,..., T and Zq — ko, Z\ — k \,..., Zt — kT are 
the risk compensating portfolio holdings after the trades at times 0,1, ...,T. Equiva¬ 
lently, the portfolio injections (respectively withdrawals - if negative) (ut)t£ t, needed 
to update the risk compensating portfolio holdings and defined by ut = Z t — Z t -i + kt-i 
(with uq = Zq), satisfy 


t -1 


u t £ A™ 


X -(- ^ ^ jus k s ) kt 


s=0 
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for every time t, and 


0 G Rt(-ut + 1 ), 

for t G T\{T}. This gives the same interpretation as with the composed regulator risk 
measure discussed in section [5] The portfolio injections of the next time period t+ 1 
are random, but acceptable with respect to the one-period risk measure Rt,t+i- 

Let us now discuss how to calculate R™ ar (X). Assume is closed and condi¬ 
tionally convex for all times t, then (R™ ar )t£T can be calculated with the approach 


discussed in sections 4.1 and 4.2 by adding K f [u; t ] to the vector optimization problem 
(3.6) at each time t. If Kf = L t (Kt) for a solvency cone Kf for all times t (as in 
example 6.1), the set-valued risk measure can also be computed directly by replacing 


the ordering cone M+ = with Kt[u>t] in problem (3.6). 

As with the ‘regulator risk measure’ considered in the previous sections, the market 
extension R™ ar (X ) might not be closed, but an arbitrarily close approximation is given 
by its closed-valued variant 


R r t nar (X) := cl J (R t (-Z) + K t ). 


zeRr+T (x) 

Solving at each node backwards in time the vector optimization problem with ob¬ 
jective 4> t (Z, Y + k ) and feasible region ( Z,Y,k ) G Z t x K t [u)t] (with Rt+i replaced 
by R™ f) yields (R] nar )teT- Benson’s algorithm yields the solution set X™ ar [u)t\ C 
L t+1 [uJt\ x R d x K t[ut] (with L t+l [uit\ = L°(U succ(cut), 2 succ ^ ut \ P(-|cjt);M d )), which can 
be used to calculate the sequence {Z t )t^T\{ o} and now additionally the sequence of 
trades (kt)te T forwards in time on the realizing path for a given Zq G R,Q iar (X ) (or 


its inner approximation) at time t = 0. Utilizing (5.4), we can find a convex com¬ 
bination of the Y elements of the solution to describe any portfolio on the efficient 
frontier of the risk measure, respectively of the inner approximation of it in the con¬ 
vex, non-polyhedral case. The same convex coefficients are then used for both the 
trading strategy and the next time step full capital requirements, which are then used 
as the starting value in the next period. 

Remark 6.4. Let us comment on the choice of M = M. d in this subsection. After 
calculating the risk measure (i?™ ar ) tg t, it is of course possible to choose a subspace 
of eligible portfolios M and choose the capital requirements to be in that space (if 
Rt iar (X) n M t is non-empty). 

However, if the subspace M ^ were to be chosen first and used for the recur¬ 
sive computation, then the last equality in (6.2) would no longer hold. This would 


cause several problems as the market extended multi-portfolio time consistent version, 


i.e. 


R 


mar.M (jy) ._ (J A , gKf UzeR™“ r (X-fc) ^t{-Z) PI M t , while retaining the capital in¬ 
jection interpretation given above, is in general not multi-portfolio time consistent. 
Furthermore, it would be more restrictive than the approach proposed above as it 
holds R™ ar ' M (X) C R? ar {X) n M t . 

On the other hand, the multi-portfolio time consistent version of the market ex¬ 
tension, while being multi-portfolio time consistent, does not admit a good economic 
interpretation for M ^M. d . 

All of these problems disappear if M = M. d . For this, and the motivation given at 
the beginning of this subsection, we suggest to use M = when considering market 
extensions and if needed one can choose ut G R 1 [ Lar {X) n Mt afterwards. 
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Example 6.5. We can use the algorithm from section [3| to compute the set of super¬ 
hedging prices under either a conical or convex market model. Dual representations 
for the set of superhedging portfolios are considered in e.g. n3[g|T30/ under pro¬ 
portional transaction costs; and in e.g. JM W under convex transaction costs. We 
calculate the set of superhedging portfolios by computing the market-compatible version 
of the worst case risk measure. Let {Rt)tei be the worst case risk measure, that is 
Rt : L —»• V[L t \ L t + ) with 


Rt(X) = {ueL t :X + u& L + }. 


The worst case risk measure (Rt)te t conditionally convex and polyhedral with Rt,t+i 


in (4.1) given as the upper image of a linear vector optimization problem. By proposi¬ 
tion f.l one can calculate (Rt)te T uit-wise since the worst case risk measure is multi¬ 


portfolio time consistent. 

Consider the market extension of the worst case risk measure, where trading is 
modeled by a sequence of solvency regions (Kt)t& t- The multi-portfolio time consistent 
market extension (Rf iar ) t £t with = L t {K t ) is nothing else than the superhedging 
risk measure. In particular, for a given claim X E L, the set SHPt(X) := Rf nar (—X) 
is the set of superhedging portfolios of X. 

Under proportional transaction costs, modeled by a sequence of solvency cones 
(K t ) te t, by proposition f.l and the discussion in section [h| the set of superhedging 
portfolios SHP t {X) of X can be calculated backwards in time by solving a sequence of 

This backwards 
which could be 


linear vector optimization problems (3.6) with ordering cone K t [oJt\. 
recursive algorithm is exactly the one proposed in m, see also J2S 
obtained as the simple structure of the worst case risk measure (Rt)te T yields a great 

Note that this 


simplification to the recursive structure ( 2 . 2 ), respectively 


(2.4). 


simplification is specific to that example, which means that the method in JWl ‘23j can¬ 
not be generalized to other risk measures, whereas the approach discussed in this paper 
is widely applicable. 


7 Numerical examples 

In this section we will apply the algorithms for the recursive calculation of polyhedral 
and conditionally convex risk measures presented in section [3j Specifically, we will 
consider the superhedging risk measure, relaxed worst case risk measure, average value 
at risk, and entropic risk measure. 

We consider a multi-dimensional tree that approximates the d — 1 asset prices 
(denoted in the domestic currency) and assume stock price dynamics under the physical 
measure P are given by correlated geometric Brownian motions: 

dS\ = SKfiidt + a t dW[), i = l,...,d—l 

for Brownian motions W l and W 3 with correlation p t j E [—1,1]. To create a tree 
for the correlated risky assets, we follow the approach in [IS ]. We expand the tree 
structure produced in such a method by allowing for n d ~ l (recombining) branches for 
any natural number n > 2 and consider some maximum change from a parent to child 
node given by v E M++ (instead of 2 d 1 branches and v = 1 from the binomial model 
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presented in [IS]). That is, since every asset can rise or fall, we consider the set of 
possible up-down scenarios given by 


£ = { (wi ,..., Wd-i) : Wi G 


2v 

-y,-, 

n — 1 


2u 

V -, V 

n — 1 


Vi = 1, ...,d — 1 


We note that in the situation with only a single risky asset, n = 2, and v = 1, this 
reduces to the Cox-Ross-Rubinstein binomial tree model. To calculate the (conditional) 
probabilities of reaching a successor node, we partition the space 1 into boxes 
in such a way that each element of £ resides in a unique box, then the probability of 
rising or falling by level e G £ is given by the probability of the (multivariate) normal 
distribution over the box containing e. 

For simplicity, we additionally assume that the proportional transaction costs are 
constant for each of the risky assets, given by 7 = (71, ...,7d-i) T G M^ 1 (possibly 0 ). 
Thus the bid and ask prices are given by (S ) 6 )* = 5|(1 — 7 *) and (Sf) 1 = ££(1 + 7 *) 
respectively for every i = 1, ...,d — 1. In the case that 7 * = 0 then the bid-ask spread 
is 0 for the i th risky asset. 

Assume the existence of a risk-free asset with dynamics (Bt)t£T and no bid-ask 
spread, i.e. B% = Bf at all times t. Further, we consider the case where cash (i.e. the 
risk-free asset) is an intermediary for all transactions. That is, the exchange between 
any two assets is done via cash and not directly. Under proportional transaction 
costs, the above simplifying assumptions ensure that the solvency cone Kt at time t is 
generated by the columns of the matrix 


B t 


-m 


-1 


(d-l)x(d-l) £(d—l)x(d—l) 


where Fd-i)x(d-i) denotes the identity matrix with d — 1 rows and columns. 

For the examples that contain convex transaction costs, we restrict ourselves to the 
two asset case (a single risky and a single risk-free asset; assets 1 and 0 respectively) 
and generate the solvency region Kt at time t by the dual equations: k G Kt if 


K+(k) := 


k 0 + 0t,c 


it, 1 


1 -exp r^ 


1 — exp ( — 


S?8 t 


S*k A 
0 

+ h 


> 0 . 


In the above equations we specify a parameter 0t G L t (M?_ + ) which defines the max¬ 
imum number of risky asset that can be bought with the riskless asset 9t,i and the 
maximum number of the riskless asset that can be bought with the risky asset Of ,0 at 
time t. That is, beginning from the 0 portfolio, it would not be possible to attain more 
than 9t ,0 units of the risk-free asset or 6t, 1 units of the risky asset. For simplicity, we 
assume that the trading strategy chosen does not impact the future market. Thus, a 
trade at time t does not affect the market at time t + 1 . 


Example 7.1. Consider a market with proportional transaction costs and two assets 
(a risk-free bond and a risky asset). We choose as a market model a recombining tree 
with 25 branches and T = 9 time steps over a one year time horizon. Further consider 
the maximal possible rise or fall in the Brownian motion to be given by v = 2. Consider 
the market with high proportional transaction costs, defined by 7 = 30%. 


19 









Let the risk-free rate of return be 10%. Let the drift for the risky asset be /a = 12.5% 
and the volatility given by a = 0.5. Consider the initial value of the risky asset to be 
Sq = $100 (measured in the risk-free asset). 

Consider the superhedging risk measure and the average value at risk with constant 
parameter A = (30%, 30%) T on the terminal payoff X of an at the money European 
put option, i.e. with strike price $100. Running the polyhedral algorithm presented in 
this paper, the efficient frontier of the time 0 superhedging risk measure and composed 

-—v—- A ,mar 

market extended average value at risk (X) is given by figure\f^ The circles 

in figure [I] denote the vertices of the risk measures. It is clear that the superhedging 
risk measure provides a more conservative set of risk compensating portfolios. Note 
that the deviation from a line to the efficient frontier of the superhedging risk measure 
is due entirely to the transaction costs, which is in contrast to the average value of risk 
as the corners there are not solely determined by the transaction costs. 



Figure 1: Example |7.1[ The superhedging risk measure and composed average value at risk 
under high proportional transaction costs 


Example 7.2. Consider a market with proportional transaction costs and three assets 
(risk-free bond and two correlated risky assets). We will approximate the market with 
a binomial tree model with T = 20 time steps over a one year time horizon. Consider 
a market with proportional transaction costs defined by 7 = 5%. 

Let the risk-free rate of return be 10%. Let the drift for the risky assets be given by 
pi = 15% and p -2 = 30%. Let the volatility for the risky assets be given by (T\ = 0.5 
and <72 = 1. Let the correlation be given by p = 0.5. Consider the initial value of the 
risky assets to be Sq = ($ 1 , $ 1 ) T (measured in the risk-free asset). 

Let X be the terminal payoff of an outperformance option with strike price K = 

$1.10, i.e. X= ^-X/{m ax (5“)>iC}^{(S“) 1 >(S'“) 2 ,(5“) 1 >K}) -^{(S“) 2 >(5“) 1 ,(5f,) 2 >X}) • 
Consider the relaxed worst case risk measure with constant parameters £i = .25 
for i = 0 , 1,2 and G is the convex cone generated by the vectors (1, —.25, —.25) T , 
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(—.25,1, —.25) T , and (—.25, —.25,1) T . The market extended multi-portfolio time con¬ 
sistent version of the relaxed worst case risk measure can he calculated via the polyhedral 
algorithm presented in this paper; a contour plot of the efficient frontier of the time 0 
risk measure Rff ar {X) is given by figure^ Notice that, as the capital in the risk-less 
asset increases the level of capital necessary in the risky assets decreases. However, 
increasing capital in one risky asset cannot totally offset a decrease in capital in the 
other risky asset, as evidenced by the curvature of the contour lines. 



Figure 2: Example |7.2| Contour plot of the efficient frontier for the composed relaxed worst 
case risk measure with 3 assets at differing levels of capital in the risk-less asset 


Example 7.3. Consider a market with convex transaction costs and two assets (risk¬ 
free bond and a risky asset). We consider the 2 time step Cox-Ross-Rubinstein model. 
We consider a market with proportional transaction costs given by 7 = 5% and convex 
transaction costs given by 9t = (500, 500) T almost surely for each time point. 

Let the risk-free rate of return be 0%. Let the drift for the risky asset be /a = 12.5% 
and the volatility given by a = 0.5. Consider the initial value of the risky asset to be 
So = $1 (measured in the risk-free asset). 

Let X be the terminal payoff from an out of the money binary option paying out 
$10 with strike price $1.20. We are able to compute an approximation of the efficient 
frontier at time 0 for the superhedging and entropic risk measures by running the convex 
algorithm presented in this text. We consider two cases for the entropic risk measure, 
each with constant parameters X\ = 10% for each asset and time and Ct = L t (C) 
where: 

1. C = cone((l,0) T , (0,1) T ) be the convex cone generated by the vectors (1,0) T and 
( 0 , 1 ) T , i.e. the restrictive entropic risk measure; and 

2. C = cone((l, -0.90) T , (-0.90,1) T ). 

Running the convex algorithm presented in this paper, with the approximation error 
at time t = 0 is given by e < 0.30 in all cases, the efficient frontier of the time 0 
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composed market extended risk measures is shown in figures 00 and 1 5| The fivst 
plot, figure 0 shows that at a large enough scale the different risk measures all appear 
identical. The curvature of the risk measures is also very evident at this size, the 
asymptotic behavior at —1500 in each asset is due to the choice of 9. In figure 0 we 
can see discrepancies between the superhedging risk measure and the two (composed) 
entropic risk measures. However both (composed) entropic risk measures still appear to 
coincide. In the final plot, figure [5| the distinction between the least restrictive entropic 
risk measure is pronounced at this zoomed-in level of detail. The superhedging risk 
measure provides the most conservative estimate, in line with the theory. 



Figure 3: Example |7.3| Convex superhedging and entropic risk measures under proportional 
and convex transaction costs, zoomed out view 
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Figure 4: Example 7.3 Convex superhedging price and entropic risk measures under pro¬ 


portional and convex transaction costs, mid-sized view 



Figure 5: Example 7.3 Convex super hedging price and entropic risk measures under pro¬ 


portional and convex transaction costs, near 0 
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